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1 Introduction 

Current sources within the brain generate weak 
magnetic fields outside the head [1], These magnetic 
fields depend on the location and the features of the 
sources. The purpose of the present study is to 
investigate the effects of extended sources located in 
deep brain structures, such as the midbrain or 
brainstem, by analytical model calculations. 

For the sources two line source configurations of 
different properties are considered. The line sources 
are assumed to be located deep in a conducting 
sphere of constant conductivity used to model the 
influence of the head. 

The magnetic field outside the spherical head model 
is studied by means of multipole expansion. This 
expansion is represented in terms of magnetic 
multipole moments at the center of the conducting 
sphere. The information content of the magnetic 
multipole fields is analyzed up to the octupole term 
for locations of the line sources near the center of 
the head model. As a special case a location of the 
source models at the center of the conducting sphere 
is considered. For this situation the magnetic field of 
a current dipole source is known to vanish [1, 2], 
The results allow general conclusions concerning 
the magnetic field of sources at the center of a 
conducting sphere. 

In addition to the forward problem, multipole 
expansion is also used for inverse calculations. For 
both source models the single current dipole solution 
of the inverse problem is derived analytically. 

2 Methods 

The magnetic field outside a conducting sphere is 
represented as the sum of the source term and the 
contribution of the volume conductor [1, 2, 3], For 
the field of the source described by a primary 
current density a spherical harmonic multipole 
expansion is performed. The origin of the coordinate 
system is taken as the center of the conducting 
sphere. The expansion of the source term is charac¬ 
terized by two types of multipole coefficients which 
can be classified as the electric moments A nm ,B nm 
and the magnetic moments oi nm ,P mn [4, 5, 6, 7], The 


magnetic field contributions of the electric multipole 
moments and the conducting sphere cancel each 
other. Therefore, the magnetic field is given by the 
expansion 
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with the magnetic multipole moments 
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In these expressions, B is the magnetic induction, 
and is the permeability of free space. Re denotes 
the real part, i is the imaginary unit, V is the del 
operator, and r,B,(p are the spherical coordinates. 
Further, Y"' (d,(p) = P"‘ (cos i})e im,p are the spherical 
harmonics, where P” (cos d) are the associated 
Legendre functions of the first kind. The bar over 
Y"'(B,(p) denotes the complex conjugate. J p 
represents the primary current density. In addition, 
§nm = (2 - S m0 )(n - m)\/(n + m )\, where 8 m0 is the 
Kronecker delta which is unity for m = 0 and zero 
otherwise. The index n gives the order of the 
multipole. Note that n = 1 is the magnetic dipole, 
n = 2 the magnetic quadrupole, n = 3 the magnetic 
octupole, etc. 

The effects of subcortical sources on the magnetic 
field outside the conducting sphere are studied for 
two line source types. Model LSI is a line carrying a 
constant current / pointing in x direction of the 
Cartesian coordinates. The spatial extent of this line 
current is chosen to be 2a. Line source type LS2 of 
extent 2b has a uniform distribution of perpendicular 
current dipoles pointing in x direction and is parallel 
to they axis. The geometrical center of both sources 
is assumed to be located at z 0 on the z axis. These 
source types represent the continuous versions of the 
in-line source and the side-by-side line source, 
respectively, of the numerical study [8] on EEG and 
MEG dipole inverse solutions for various cortical 
sources in a spherical head model. 



In addition to the magnetic field calculations, 
multipole expansion is also used to determine the 
single current dipole solution of the biomagnetic 
inverse problem. To obtain this equivalent current 
dipole (ECD) a method developed for the electric 
potential [9, 10] is applied to the magnetic field. The 
ECD is considered to be a source which generates 
the same magnetic dipole and quadrupole fields of 
the multipole expansion like the actual line sources 
located near the center of the conducting sphere. The 
higher-order terms beyond the magnetic quadrupole 
are neglected. The ECD parameters are obtained 
from a system of equations resulting from the 
magnetic dipole and quadrupole moments of the 
multipole expansion. In the present case this system 
reduces to the equations for an ECD of the current 
dipole moment p = ( A U ,0,0) located at z on the z 
axis. They are given by 
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with the solutions 
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where /?,, and /) 21 are the magnetic dipole and 
quadrupole moments, respectively. 


3 Results 

3.1 Magnetic field 

For the magnetic field calculation on the basis of 
Eq. 1 it is necessary to know the magnetic multipole 
moments which are represented by Eq. 2. 

In the following for both line source types the 
resulting magnetic moments are presented up to the 
octupole term. The geometrical center of the sources 
is assumed to be located at z 0 on the z axis outside 
the center of the conducting sphere. 

First the line source LSI is considered. In the order 
n = 1 Eq. 2 gives the magnetic dipole moment 


For n = 3 the resulting magnetic octupole moments 
are given by 
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and 
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In the case of the line source LS2 the magnetic 
dipole moment is 
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where A u =2 qb is the electric dipole moment of 
the line source of constant current dipole density q 
and extent 2b. The magnetic quadrupole contribu¬ 
tion is given by 
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The magnetic octupole moments are found to be 
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For both line source types the magnetic multipole 
moments up to the octupole contain all information 
to determine the location of the geometrical center 
and the extent of the sources. In the limits a —> 0 
and b —> 0, respectively, these multipole moments 
reduce to the magnetic multipole moments for a 
single current dipole of the electric dipole moment 
A u located at z 0 on the z axis [11]. 

The magnetic field calculated by means of the above 
magnetic multipole moments has for both line 
source types the form 


Pn = 


z 0 A n 


(5) 


where A u = 2aI is the electric dipole moment of 
the line current I of extent 2a. The next contribution 
which is obtained for n = 2 is the magnetic quadru¬ 
pole moment 


0 #* 


z';A.. 


( 6 ) 


B(r) = B n=l (?) + B n=2 (?) + S„ =3 (?) + •••. (13) 

Using Cartesian coordinates to present the contri¬ 
butions to the magnetic field the magnetic dipole 
term is given by 

B n=l (4kA0uOyr-r 2 e y ), (14) 
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where e x ,e ,e z are the unit vectors. The magnetic 
quadrupole term is found to be 
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Discussion 
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The contribution of the magnetic octupole can be 
shown to be 

B n= ,{r) = ^\W iX [(lz 2 -r 2 )5yr 
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Now it is assumed that the geometrical center of the 
line sources is located at the center of the conducting 
sphere, i. e., z 0 = 0. For the line source LSI the 
magnetic multipole moments up to the octupole 
presented above and all higher-order magnetic 
moments vanish. Therefore, for the line current LS1 
the magnetic field B is equal to zero. In contrast, 
for the line source LS2 the magnetic field does not 
vanish. In this case the magnetic quadrupole term is 
found to be the lowest-order contribution to the 
magnetic field. 

3,2 Current dipole inverse solution 

The location and the moment of the ECD for the 
magnetic field of the line sources located near the 
center of the conducting sphere are obtained from 
Eq. 4 and the expressions for /f M and /J 21 presented 
above. 

In the case of the line source LSI the ECD result is 
given by 

z = z 0 ,A n =A n . (17) 

Therefore, for the line current the ECD represents 
the location of the geometrical center on the z axis 
and the electric dipole moment. 

For the line source LS2 the ECD solution is 

z = kz 0 , A n = (18) 

k 

with 

k = \ ——-T-. (19) 
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This result describes the effect of the source extent 
on the location and the moment of the ECD. In 
contrast to the line current LSI, for the source LS2 
the inverse solution shows systematic deviations 
from the actual source parameters. 


The present study provides insight into the magnetic 
field from extended sources at deep locations in a 
spherical head model. The analytical results ob¬ 
tained by means of multipole expansion written in 
terms of magnetic multipole moments at the center 
of the conducting sphere allow a detailed analysis of 
the effects of the sources. 

First the source center is chosen to be located at z 0 
on the z axis outside the center of the conducting 
sphere. The results show that for the line current 
LSI the first two multipole terms up to the magnetic 
quadrupole do not depend on the source extent. In 
contrast, the magnetic quadrupole term of the line 
source LS2 depends on the source extent. For both 
line source types the octupole contributions are 
influenced by the source extent. Especially, the /f 33 - 
term is found to be caused by the extent of the line 
sources. This term vanishes in both cases with 
vanishing source extent. 

Now the source center is taken at the center of the 
conducting sphere, i. e., z 0 = 0 . For the line source 
LSI the magnetic field vanishes, because all 
magnetic multipole moments are equal to zero. 
However, for the line source LS2 the magnetic 
induction is given by 

B(r)= 3 /° 7 P ?l (5yzr - zr% -yr 2 e z ) + - (20) 
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with the magnetic quadrupole moment 



Thus, the magnetic field does not vanish for the line 
source LS2 with its geometrical center at the center 
of the conducting sphere. This significant result is 
valid for all source types for which the multipole 
source description gives nonzero magnetic multipole 
moments [12]. 

Furthermore, in this study the multipole expansion 
of the magnetic field provides the basis for inverse 
calculations. In analogy to the electric case [9, 10] 
the magnetic dipole and quadrupole moments are 
used to determine the ECD for deep locations of the 
actual sources near the center of the conducting 
sphere. The magnetic octupole contribution and all 
higher-order term s are neglected. In the case of the 
line source LSI the ECD describes the location of 
the center and the electric dipole moment of the 
actual source. The reason is that the first two 
multipole terms do not depend on the source extent. 
For the line source LS2 the ECD result shows 



systematic deviations from the source parameters z 0 
and A n , because the magnetic quadrupole moment 
is affected by the extent of the actual source. 

For the line source LS2 the location of the source 
center z 0 on the z axis outside the center of the 
conducting sphere and the electric dipole moment 
A n can be determined by taking into account the 
next higher-order term of the multipole expansion, 

i. e., the magnetic octupole. The results for the 
magnetic multipole moments up to the octupole 
contributions presented above can be used to derive 
the equation 

2/J n Zp - 3 jS 21 z 0 +32^3 = 0 . (22) 


The location of the source center z 0 on the z axis is 
given by the solution 
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This result can be used to determine the electric 
dipole moment A n of the line source by means of 
the relation 



In addition, the extent of both line sources can be 
obtained by taking into account the magnetic 
octupole term. 

For the determination of the line source parameters 
it is necessary to know the corresponding magnetic 
multipole moments which can be calculated from 
the magnetic field outside the conducting sphere by 
using a surface integration method. Considering the 
radial component of the magnetic induction the 
magnetic multipole moments can be obtained by 
means of the equation [7] 

a nm + ifinm = 2n + lgnm f r n Y n m (i &,<p)B(r)dS , (25) 
n +1 H 0 i 

where S is a spherical surface enclosing the con¬ 
ducting sphere. 
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